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^ ■ I. INTRODUCTION 

' The muon decay lifetime (r^) has been used for long as an input parameter for high precision predictions of the 
Standard Model (SM). It allows for an indirect determination of the mass of the W boson {Mw), which suffers 
currently from a large experimental error of 39 MeV [ ll , one order of magnitude worse than that of the Z boson mass 
I , (Mz)- A reduction of this error by LHC to 15 MeV pj and by a future linear collider to 6 MeV |^ would provide a 
^ • stringent test of the SM by confronting the theoretical prediction with the experimental value. 

[ The extraction of Mw with an accuracy matching that of next experiments, i.e. at the level of a few MeV 
00 ' necessitates radiative corrections beyond one loop order. Large two-loop contributions from fermionic loops have 

■ been calculated in The current prediction is affected by two types of uncertainties. First, apart from the still 
unknown Higgs boson mass, two input parameters introduce large errors. The current knowledge of the top quark 
mass results in an error of about 30 MeV [||, which should be reduced by LHC to 10 MeV and by a linear collider 
even down to 1.2 MeV. The inaccuracy of the knowledge of the running of the fine structure constant up to the Mz 
scale, Aa{Mz), introduces a further 6.5 MeV error. Second, several higher order corrections are unknown. In fact 
the last unknown correction at the C(a^) order has been calculated only recently in and 0. This contribution 
comes from diagrams with no closed fermion loops. 

O I It is the purpose of the present work to give a detailed description of the methods used in the calculations presented 
1^ in 1^ and [Q. Since one of the groups used high precision numeric methods and the other deep expansions both in 

mass differences and in large masses, a comparison can be given. 
J> , In the next section we discuss the question of matching of the Fermi theory onto the Standard Model at low 
energy scales. Then we move to the discussion of renormalisation in the on-shell scheme and continue with the MS 
^ . , scheme. A section on the transition between the schemes contains comparisons of the methods and the final results. 

■ The description of computational methods and conclusions close the main part of the work. In the appendices, a 
" ■ ■ derivation of the electric charge counterterm through the U{1) Ward identity can found, followed by the explicit 

analytic results of the expansions of the on-shell and the MS quantities. 



o 

O 

Oh: 



X 



II. MATCHING 



The muon lifetime can be computed from the effective Fermi theory given by the lagrangian 

£cff _ Op + higher dimension operators , (1) 

v2 

where Op the 4-fermion Fermi operator of dimension six 

Of = [P^ni - 75)/^] [e7a(l - l5>e] , (2) 

and is the Fermi constant. Note that Eq. ^ is a definition of Gp. This lagrangian can be used to describe low 
energy processes (such that energies are <C Mw) mediated by the weak charged current. Since the theory Eq. is 
nonrcnormalisable, an ultraviolet cut-off A should be introduced. 
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In particular for the muon decay process we have 

1 _ G|.m5 



1927r^ V ml 



(1 + Ag) , (3) 



with TOe and to^ being the masses of the electron and the muon respectively. The quantity Aq describes all QED 
corrections in the Fermi theory and has been calculated at the one-loop Q and at the two-loop ||^ order. 

By its nature G_f is the Wilson coefficient function of the operator Of and can be evaluated from the SM. Tradi- 
tionally the matching relation between Gp and the parameters of the SM is parametrised as follows 

Gp 

The quantity Ar — Ar^^^ + Ar*^^) + . . . absorbs the effects of all loop diagrams. 

It is the purpose of the following subsection to establish the framework for the calculation of Ar. 



A. Factorisation theorem 



In principle, the muon decay amplitude can be evaluated directly in the SM, but this is not feasible in practice. 
There are many scales involved which vary from less than 1 MeV to 100 GeV, i.e. by more than 5 orders of magnitude! 
On the other hand the number of Feynman diagrams grows very fast with the number of loops. A way out to keep 
the problem manageable is to switch on the machinery of effective lagrangians (see Eq. This allows one to simplify 
the calculation enormously and to separate consistently the low energy ( "soft" ) dynamics from high energy ( "hard" ) 
static characteristics. 

Suppose that we can compute the muon decay amplitude ^j-^g gjyj q;^]2en the Fermi constant Gp defined 

through Eq. |l| can be predicted from A^^. Indeed we should require that both evaluations in the SM and the Fermi 
theory give the same result. At the tree level the corresponding matching equation reads 

A- = e|(,|0He.„..) + 0(^). (5) 

This equation just states that the amplitude of the process /i —f evv is the same both in the full SM and in the 
effective Fermi theory up to operators of higher dimension. 

When loop effects are taken into account, matrix elements in both sides of Eq. ^ get quantum corrections. Since 
and {[i\Op\ev^J7^ are amputated matrix elements one has to renormalise also the external wave functions. Therefore 
the final form of the matching equation reads 



cff 7cff 7-1 '^P .T \ I f ™^ 



where and Z'^^ are wave function rcnormalisation constants of the fcrmions evaluated in the SM and in the 
effective theory respectively and Zof. is the rcnormalisation constant of the Fermi operator in the effective theory. 
There are two ways to compute Gp from the SM: 

1. standard matching calculation, or 

2. automatic matching via factorisation theorem. 

The former approach works always by simply computing all ingredients (apart from G^) in the matching equation 
Eq. 1^. This requires however much extra efforts to evaluate the "soft" pieces (or, at least, to separate them) in 
the amplitudes and the Z's. Historically for this purpose the Pauli-Villars regularisation was used in [|o) and then 
extended to two- loop order in . The same approach has been applied also in js] . 

How it works at the 1-loop level is demonstrated in Fig. |l|. There are only three infrared divergent diagrams with 
photon. From each diagram its counterpart in the Fermi theory should be subtracted. The left diagram in each line of 
Fig. corresponds to the result in the full model and therefore contains both the "soft" and the "hard" part. The right 
one contains only the "soft" part, which means that the difference is the requested "hard" correction. In addition, for 
the diagrams in the frame the Pauli-Villars regularisation is introduced to regularise the ultraviolet divergences. At 
the two- loop level we have a very similar situation. The difference is that instead of "hard" and "soft" terms there 
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FIG. 1: one-loop factorisation with Pauli-Villars regularisation. 

are now "hard-hard", "hard-soft", "soft-hard" and "soft-soft" contributions. From these only the "hard-hard" piece 
contributes to Gp- 

Accidentally, it happens that the sum of the three "soft" diagrams inside the frame in Fig. |^ is an ultraviolet finite 
quantity (let us call it Sgoft)- It is easy to prove that this holds true also to all orders. This is a consequence of the 
Ward-Takahashi identity for QED. This fact, however, is a pure coincidence rather than something fundamental. If 
such a cancelation had not occurred, renormalisation of the operator Op would be required as it is taken into account 
in Eq. ^. 

The scheme given in Fig. |^ is consistent but the disadvantage of it is that there arises the problem of bookkeeping 
of "soft" and "hard" parts and already at the two-loop level the problem becomes very complicated. Indeed, at the 
two- loop level one has to subtract from each diagram the "hard-soft", "soft-hard" and "soft-soft" pieces. 

Therefore it would be very helpful to find some other way to obtain the "hard" part. Thus we come to the second 
way to compute Gp — automatic matching. This procedure is the most straightforward and the most economical 
(minimal in costs) way to compute. It is based on the factorisation theorem, proven e.g. in It allows one to 

extract the "hard" part directly without any reference to "soft" pieces. As a well known example of such a procedure 
we can mention the evaluation of Wilson coefficient functions in deep inelastic scattering processes. 

Returning to the sum of the three "soft" graphs in Fig. |l] (Sgoft) we notice that in Gp all "soft" modes are eliminated. 
This means, that all subgraphs in Fig. |l] should be computed at vanishing masses of the leptons. In this case the 
Ward-Takahashi identity not only makes Ssoft ultraviolet finite but also nullifies it. Thus all "soft" parts add up to 
zero. This is also true to all orders of perturbation theory. In other words, one can from the very beginning nullify 
all external momenta and masses and evaluate the obtained bubble diagrams. Of course new infrared divergences 
are generated. They cancel however in the expression for Gp. To regularise these infrared divergences we use the 
dimensional regularisation. 

To prove rigorously that infrared singularities indeed drop out from the result one can turn to the framework for 
construction of effective low energy lagrangians given in [|l2| . At the level of individual Feynman diagrams one can 
separate "soft" and "hard" scales with the help of the asymptotic expansion procedure |Q . Let F denote a Feynman 
diagram. Then 

HCP 

where the sum runs over all "hard" subgraphs H of the diagram F; 5 is a "soft" subgraph obtained from F by 
shrinking iJ to a point and T stands for the Taylor expansion (before integration!) of H with respect to all "soft" 
parameters. The exact rules for construction of hard subgraphs are discussed in details in p3| . 

The important property of the operation Eq. ^ is that it has the combinatorial structure of the _R-operation . 
This allows one to promote the operation on a single Feynman diagram to the operation on the whole Feynman 
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amplitude (the factorisation theorem). By this procedure all infrared divergencies are absorbed either by the "soft" 
matrix element or by the renormalisation constant Zq of the operator. The detailed discussion can be found in |p^ . 
In the case of Gp we have further simplifications. 

• The anomalous dimension of the Fermi operator Op \b zero, therefore Zq^ in the matching equation Eq. ^ is 
equal to one. 

• At zero lepton masses and external momenta all Z^ and the "soft" matrix element in Eq. ^ are equal to one. 
Finally we get 



V2 



7SM 7SM 7SM 7SM 



SM 



(8) 



hard 



where the subscript "hard" means that all "soft" scales are put to zero. 

Thus the problem is reduced completely to the vacuum Feynman diagrams of one- and two-loop order and the 
bookkeeping problem does not arise at all. The wave function renormalisation constants are to be computed in the 
on-shell scheme. Again, for massless leptons, the wave function renormalisation constants are defined through vacuum 
diagrams only. Such diagrams can be evaluated analytically using reduction formulae of based on integration by 
parts identities 



B. Projection 

An important problem in the calculation is the reduction of the amplitudes to scalar integrals. It is not only of 
practical importance. In fact it is connected to the correct definition of the matrix elements in the model, since 
dimensional regularisation is used. 

The matching onto the Fermi theory with its double V ~ A cliiral structure is made possible because of the left- 
handedness of the charged current in the Standard Model. The "hard" components of the diagrams contain only 
massless fermions and therefore formally the structure of the two spinor lines can be mapped onto the operator 

7^-Pl®7p-Pl- (9) 

In four dimensions, every string of an odd number of gamma matrices and a left-handed projector can be reduced to 
the structure due to the Chisholm identity 

i^iiuip = gt,vip + 9upip - gppiu - i^t^vp^i^iz- (10) 

The reduction leads to the operator 

T^.7''Pl®7"^l, (11) 

where T^^ is some tensor made of the integration momenta. Since there are no non- vanishing external momenta, this 
tensor must be proportional to g^i, and the result Eq. ^ follows. A suitable way to obtain directly the right value is 
to use a projector made of trace operators. Let the original product of strings of gamma matrices be denoted by 

Ti®T2. (12) 

We wish to obtain the proportionality coefficient A in the following equation 

Jri®T2^Ax{^''PL®J^PL). (13) 

Two possibilities of closing the spinor strings with trace operators are depicted in Fig. ||. The left one has been 
used in Q and is given by the equation 

^ = ^ / Tr(ri7'^PB)Tr(r27M^'fl), (14) 

where the dimension of space-time d has been kept arbitrary and the trace of the unit matrix has been put to 4, as 
usual. A second possibility which we used to perform the calculations presented in this work is given by 



A = -^^/Tr(r,7.P.r.7^P.), 



(15) 
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FIG. 2: Two of the possible projectors for Ar. The dashed hnes represent the strings of Dirac matrices, while the crosses, the 
projection operators. 

and corresponds to the right picture in Fig. |^ 

Both projectors are obviously equivalent in four dimensions due to the Chisholm identity as explained above. 
The difference starts to be important for divergent integrals. In fact the problem does only occur for one-particle- 
irreducible four-point diagrams, where the divergence can come from two sources. First from the external wave 
function renormalisation, which is incoinplcte due to infrared divergences and second due to infrared divergences of 
the diagrams themselves. As noticed in Q the first projector Eq. |l^ needs to be corrected, as it does not fulfil several 
requirements, like for example the vanishing of diagrams with propagator insertion in the photon lines. Moreover, one 
can explicitely check that without corrections the subtracted diagrams in the Pauli-Villars approach do not cancel and 
the dependence on the A scale remains. In the automatic factorisation approach this shows up through an incomplete 
cancellation of divergences. Notice, however that the result is gauge independent, thus it is only the finitencss of the 
result that shows that the projector is incorrect. 

On the contrary the projector Eq. [l^ does not require any corrections. It does fulfil all of the algebraic requirements 
and also yields a finite result as well as the exact cancellation of the subtraction diagrams of Fig. |] in d-dimensions 
and in all orders of perturbation theory. This useful property follows from the fact that this projector respects the 
Fierz symmetry in d-dimensions. One can check explicitely that for example 

If^l.lpPL ® Yl'-f^'PL ~ 7M7^7p7"7''-Pi ® YPl, (16) 
where means equality after projection. 

III. ON-SHELL RENORMALISATION 

Two-loop calculations within the on-shell renormalisation scheme require the knowledge of several counterterms. 
At the very least charge and mass counterterms are needed. In this section we first discuss the problem of gauge 
invariance in connection with tadpole diagrams. We then give specific expressions for the required counterterms. 

A. Tadpoles and gauge invariance of counterterms 

It has been known for a long time that the inclusion of tadpoles is necessary to obtain gauge invariant counterterms. 
In fact this property has been first noticed ||T^ shortly after the proof of renormalisability of gauge theories. A general 
proof of the Quantum Action Principle, which has for consequence the gauge invariance of on-shell processes in the 
bare lagrangian, requires the inclusion of even those tadpoles which would be cancelled by normal ordering (one loop 
tadpoles) p8| . There are, however, two disadvantages of using tadpoles in actual calculations. First, this requires the 
inclusion of diagrams, which drop in the final result. Second, one-particle-irreducible (IPI) Green functions cannot 
contain tadpoles. As long as we wish to obtain results at the least cost and by using automated software, it is 
interesting to consider alternative possibilities. 

It turns out that it is possible to prepare the bare lagrangian in such a way, that the only gauge dependent quantities 
would be the wave function renormalisation constants and the vacuum renormalisation constant, and still all of the 
tadpoles would be cancelled. Let us start by considering a lagrangian in which the bare coupling and masses are 
defined through physical processes. The masses can be equivalently defined through the position of the poles of the 
physical ^-matrix in the complex plane as recently proved In such a case all of the bare parameters would 

be gauge invariant, because they would fulfil equations that have this same property. It is important to supply a 
condition on the vacuum expectation value of the bare Higgs field vq that would resum terms of order 0{a'^). A 



6 



choice which is still consistent with gauge invariance is 

1 (\ 



2^0 (^2^o^o-M5 j 

where Aq and /io are defined through the Higgs lagrangian 

^Higgs = - iAo($t$o)2, (18) 

and $0 is the Higgs doublet. Eq. |l^ implies the vanishing of the linear term in the lagrangian. Although this term 
will be subsequently altered, the tree level contribution will always vanish. 

We now introduce an additional renormalisation of the bare vacuum expectation value 



Wo — > Wo 



Zll\ (19) 



The renormalisation constant can be used to cancel the tadpoles recursively, which implies together with Eq. 
that the first non- vanishing term in its perturbative expansion starts at order Oipi). The linear term in the Higgs 
field can now be written as 

_ y^^o ^ _^^kl^3L(^Mlfzll''{Z, - 1)H° (20) 
eo 



where the following relations have been used 



.0 - (21) 



eo 

Mo = {M^nf, (22) 



^ I ^ '''' 



At the tree level the contribution is zero, since then zi^^ = 0, as noticed above. To one-loop order, the relation 
between the tadpole diagrams and the vacuum expectation value is simple 

SziJ^ ^ 1 ^n(l\ (24) 

where jH^-* is the sum of IPI one-loop tadpole diagrams of the Higgs field. The situation gets much more complicated 
at the two-loop level 



- . , I ,,, n^) - ( 6Z^^ + SZ^^^ + 2^ + ^ + 2'-^ - 5Zp] . (25) 

sm0wMwMjj " 2 " \ " " Mjj sm9w I 



At this level an insertion of this counterterm reproduces all of the tadpole diagrams that would be included in the 
usual approach. An example is depicted in Fig. ^. An insertion of Zy into the W boson self-energy a) , leads effectively 
through the first term in Eq. to an insertion of a one-loop tadpole with with a vertex counterterm b). This 
counterterm also contains a correction to the Vacuum expectation value of the Higgs field, which reproduces the 
tadpole diagram c). 



B. On-shell scheme counterterms 



The on-shell renormalisation scheme is defined by the requirement that the masses be identified through the poles 
of the physical S'-matrix (as the real part of the pole) , while the electric charge coincide with the value measured in 
the Thompson scattering process as for example in the quantum Hall effect. These conditions are enough to fix all 
of the free parameters of the SM with minimal Higgs sector (neglecting the CKM matrix and the strong coupling 
constant). The counterterms have been given by many authors. The peculiarity of the present work is the specific 
definition of the bare masses which are gauge invariant without including tadpole diagrams. This, however, implies 
that the formulae defining the counterterms will be slightly different. 




a) 



b) 



FIG. 3: Vacuum expectation value counterterm insertion into the W boson self-energy a), reproducing tadpole insertions b) 
and c). 



At the one-loop level, the mass counterterms are related to the on-shell self-energies through 



5M, 



2(1) 
H 

2(1) 



n 



(1) 



5Mf^ = -n^p,AMl)~Ml5Zl^\ 



(26) 
(27) 
(28) 



where iliu^T denotes the transverse part of the self-energy diagrams of the boson i. For bosonic corrections to the 
Higgs boson mass counterterm the real part has to be taken due to the possible decay into aW or Z boson pair. To 
one-loop order this still yields a gauge invariant result for the renormalised amplitude. The W and Z bosons do not 
require such a treatment neither at one nor at two-loop order. 

At the two-loop order, only W and Z boson mass counterterms are needed, and they assume the form 



SM, 



2(2) 
W 



SZi^\MlSzl^^+dA4'^), 



SM: 



(29) 
(30) 



The last term in the Z boson mass counterterm, which does not occur in the W boson mass counterterm, has its 
origin in the mixing between Z and 7. If the self-energies have imaginary parts, then suitable additional terms have to 
be included as described in [||. The above formulae are valid only if the subdivergencies in the two- loop self-energies 
are renormalised. They also require the wave function rcnormalisation constants of the bosons 



szi'^ = n 



n w w T {My^ ) , 



WW,T 

(1)' 

ZZ,T 



(M^z), 



(31) 
(32) 



and the mixing rcnormalisation 



-yZ 



(33) 



The last two constants form part of the 2x2 rcnormalisation matrix of the neutral bosons 



7I/2 1 7 
^zz 2^72 
,1/2 



(34) 



The remaining two rcnormalisation constants define the photon field and can be obtained at zero momentum transfer 
from the following formulae 



5Z^1 = - 



'Zi 



<.^,'t(0)- 



.t(0), 



(35) 
(36) 
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The electric charge counterterm can be obtained in two ways. The first consists in simply calculating the scattering 
of fermions off real photons, i.e. at zero momentum transfer. This however introduces unnecessarily three point 
functions. A second possibility is to use the U{1) Ward identity. The suitable relation between the wave function 
renormalisation constants of the photon and the Z boson has been proved in |Q using the BRS symmetry. A simpler 
proof is given in appendix The one and two-loop counterterms in the on-shell scheme are given by 



J^(i) - 1 (1) I sin^M/ ^ (1) 



(37) 



The two-loop wave function renormalisation of the photon is given by the short formula 

5Z(-^=Ii^^%^{i))~\{6Z^'y)\ (39) 
whereas in the mixing counterterm, the vacuum expectation value correction makes again its appearance 

< = -^<z{0) - \^^fz^< - -^JZ^z>^r - SZi^'SZ^P^. (40) 

In the on-shell calculation the ghost sector was also renormalised. The respective constants are as in |^ up to 
an unimportant renormalisation of the ghost wave functions, the difference being dictated by simplicity. The wave 
function renormalisation constants of the ghosts and Goldstone bosons have been left unspecified. For the ghosts, these 
constants cancel trivially within every closed loop. With the Goldstone bosons, the situation is more complicated, 
since the fact that the gauge fixing term should not be renormalised induces Goldstone wave function renormalisation 
constants in the ghost sector. These can only cancel in gauge invariant quantities. This indeed happened for all the 
mass and coupling counterterms and for the complete result. 

IV. MS RENORMALISATION 

In this section we describe in detail the renormalisation of Ar in the MS scheme. Ar is computed through the 
matching procedure described before in Section ^ Here we chose the strategy of multiplicative renormalisation. 
After multiplication by the on-shell wave function renormalisation constants of external fermion fields, the result is 
expressed in terms of bare masses and bare electric charge. In order to get the MS renormalised result for G f one 
needs to substitute all bare parameters in the form 

eo = ^^Zee(/i), 

(m0)2 = Z^^mf{^i), (41) 

where e(/i) and mi{^) are the MS charge and masses respectively and /i is the MS parameter. The MS renormalisation 
constants will be specified in the next two subsections. 

Let us stress that in Eq. ^ we renormalise only the physical parameters and no renormalisation of the unphysical 
sector (ghost sector and gauge fixing parameters) is required. The renormalisation of the boson particles' wave 
functions is also not needed since it cancels anyway in the final expression. 

A few words should also be said about tadpole diagrams, which should be added in a proper way in order to 



obtain a gauge invariant result in the SM. Unlike in the approach described in the Section HI A, where the new 
counterterm Z„ for the Higgs VEV has been introduced, here we include the tadpole diagrams explicitely. This makes 
our renormalisation constants Z's from Eq. ^ gauge invariant. 

Below we present the analytical expressions for charge and mass MS renormalisation constants, needed in order to 
obtain a finite expression for Ar in the MS scheme. 

A. Coupling and masses renormalisation 

The bare charge Cq and the MS charge e are related via 



eo 
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where the constants Z's, as we shall see in the following, can depend on sinOyy- 

There are two ways to determine the MS renormalisation constant in this expression. One is to use the Ward- 
Takahashi identity given in appendix ^ to express it in terms of the gauge boson wave function renormalisation 
constants and the renormalisation constant for sin 9w 

1 = ZJ V^+o— Tif'JZz^. (43) 



2 cos J 

Then at the one- and two-loop order for on-shell charge renormalisation constants, introduced above, we have 

<L ^ -IsZiS - + (Szi^^ + i (.Z...)^ ^ (45, 

Here SZ'-^^^^ are one- and two-loop on-shcU field renormalisation constants, expressed via bare quantities. We can 
rewrite Sz!fy and Sz'^^ in terms of self-energy diagrams 

sz^ii = (<iT(o) + ^-^zgnW ^(0) + i^zwnw ^(0) + i^^z^^ln^ ^(o)) - ^sz^'lsz'il, (46) 

where this time all of the self energies are unrenormalised. All other one- loop field renormalisation constants were 
defined before in Section III B| . At the end we have an expression for the on-shell charge renormalisation constant 
expressed via bare charge, Weinberg angle and masses. Now, rewriting the bare quantities in terms of MS ones with 
yet unknown coefficients in Eq. ^ an d requiring that transition between on-shell and MS charge should not contain 
divergencies we easily extract the MS charge renormalisation constants. 

Alternatively, the renormalisation group analysis can be applied. In order to find Z^ we differentiate Eq. (p^ ) w.r.t. 
log/i^ and take into account that 

rl r 

= -T:e + f3e, (47) 



d log fi^ 2 



where 



/^e = ^^i+(^^i + -.. (48) 

is the /3-function. Since {d/d\og fi^) Cq = 0, the l.h.s. of Eg. becomes zero after the differentiation, while the r.h.s. 
relates the coefficients bj and the unknown constants in ( [42|) 





= ^0, 


7(2,2) 


- ^6^ 


^(2,1) 





(49) 

The function (3^ can be extracted from the existing calculation in the unbroken theory. Namely, for the SU{2) and 
U{1) charges g and g' respectively, the /3-functions read 



1 g'-' , 1 g'' 3 g'^g^ 



12 167r2 4 (167r2)2 4 (167r2)2 



^ 43 g' 259 g^ 1 g^g'^ 

12 167r2 12 (167r2)2 ^ 4 (167r2)2 ' 

The one- loop result is given in while the two- loop coefficients have been evaluated in p^ . 
From the relation 

1 _ 1 1 
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FIG. 4: The exact result (solid line) for the Ar^^^ in the in the MS scheme (left and right panel) vs. its 1) large Higgs boson 
mass expansion (long dashed line) and 2) mass difference expansion (dotted line). The short dashed line 3) represents the [3/3] 
Fade approximant. In the right pannel the large Higgs boson mass expansion curve coincides completely with the numerical 
result for this range. 



it is easy to deduce that 



/?e = e^(f| + f|)- (52) 



Using now Eqs. ^ and 52 we obtain 

7_e 

'2 167r2 ' (167r2)2 \^ 6 sin 



7 e3 / 125 1 , 



and, finally, from Eq. we have 



7 

^2' 
147 



12 sin^^vK 4cos20v^- ' 



The explicit calculation confirms the above result. 

Similarly to the charge renormalisation we write for the masses of the Z, W and the Higgs bosons 

i^v) -^Vif^) ^Q^2,ZV + ^lQ^2y,Zv + (16^2)2,2 )■ (55) 

For mz and mw the renormalisation constants up to two loop are required while for the Higgs boson we need only 
the one loop expression. The analysis, similar to that described above for the charge, has been done in details in P3[ . 

There the explicit expressions for Zy'^\ Zy'^^ and Zy'"^^ are given. 



B. MS results for Ar 

In Fig. 4 we plot ArJ^^(, as a function of the MS Higgs boson mass in different scales. As input parameters we 
used the on-shell values given in Table I. 

The solid curve represents the exact result. Two other curves represent expansions in different regimes: as niH 
mz and as mn — > oo. They cover almost the whole region of the mn under consideration. In order to extend the 
range of the expansion around mz the Fade approximant has been constructed. It sufficiently improves the situation 
for the intermediate Higgs boson masses. Thus the expansions cover completely the region of interest. The details of 
expansions are discussed more precisely in Section VI B. 
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a 

Mw 
Mz 



137.03599976(50) 
80.423(39) GeV 
91.1876(21) GeV 



TABLE I: Parameter values used in the calculation jl|. 



TRANSITION BETWEEN THE SCHEMES 



Once we have the result in the MS scheme it is necessary to translate it into the on-shell parameters, which are 
known with high precision for the electroweak sector contrary to the strong interacting sector of the Standard Model. 
To this end one has to consider the proper scheme independent quantity which is 



5 1 + Ar = 

2M^ sin^ V2 



(56) 



This should be contrasted with the naive approach of taking simply Ar and substituting MS parameters. 

Using the methods described in Section IV, we obtain the following series expansions connecting on-shell and MS 
parameters 
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(58) 

(59) 
(60) 



The series for the Higgs boson mass relation is only needed to first order, since the Higgs field starts to contribute to 
the decay only at the one-loop level. 

The above relations have to be inverted to yield MS parameters in terms of the on-shcU ones. For any parameter 
A the relation will be written as follows 
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(61) 



The expansion coefficients are obtained by inverting the original series up to the required order, 
leads trivially to 



At one-loop this 
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(62) 



The coefficients for the three bosons are depicted in Fig. ^ with parameters values as given in Table |, in a comparison 
of the different evaluation methods. For Higgs boson masses greater than 200 GeV the large mass expansion with six 
coefBcients is indiscernible from the numeric result. The mass difference expansion fails always around 120 GeV. In 
the visible range from 80 GeV to 200 GeV, the Fade approximation based on the mass difference expansion turns out 
to practically coincide with the exact result for the vector bosons. For the Higgs boson this cannot happen due to 
the occurrence of the two-particle production thresholds and indeed there is a region between the thresholds which 
cannot be reproduced with neither the mass difference nor the large mass expansion. Obviously if it was needed this 
region could be covered by threshold expansions. 

The two-loop correction contains terms coming also from the one-loop terms and the proper expression reads 
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(63) 



M,,os 



The corrections for the vector bosons are depicted similarly to the one-loop case in Fig. I The expansions themselves 
are less precise. It is however interesting to note that the Fade approximation together with the large mass expansion 
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FIG. 5: One-loop corrections to the relations between the on-shell and the MS masses for the W, Z and Higgs bosons 

(AJ — a/{ATT)Xl). The long dashed lino 1) represents the largo Higgs boson mass expansion, the dotted line 2) represents the 
mass difference expansion. The short dashed line 3) gives the [4/4] Pade approximant which coincides for this range with the 
exact result for the W and Z boson mass corrections. 
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FIG. 6: Two-loop corrections to the relation between the on-shell and the MS masses for the W and Z boson (A? = 
{a/{4-iT))^Xl). The long dashed line 1) represents the large Higgs boson mass expansion, the dotted line 2) represents the 
mass difference expansion. The short dashed line 3) gives the [4/4] Pade approximant. 
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FIG. 7: On-shell Ar^J^,. The long dashed line 1) represents the large Higgs boson mass expansion, the dotted line 2) represents 
mass diflterence expansion. The short dashed line 3) gives the [3/3] Pade approximant. The dash-dotted 4) and 5) lines 
correspond to lower terms in the large Higgs boson mass expansion, whereas 6) is the leading term. 



cover the whole range with high precision. Even the threshold region is reproduced with a relatively small error, 
although this is due to the fact that the peaks are not very pronounced. 

We can now combine all the perturbative expansions and translate the MS result into the on-shell one. We shall 
not reproduce the formula since it can be easily obtained from the previous equations. It is important however to 
note two things. First, in the expression for the two-loop Ar there are the following terms 



f Ar(^)) °' = . . . + - + ... (64) 

V / sin^Ow sin^^Jw ^ ^ 



If this is combined with the fact that the results in both MS and on-shell schemes behave as 1/ sin'' 6w, it is obvious 
that one term in the W and Z boson mass difference expansion is lost. Second, the result in the MS scheme behaves 
as M'^, whereas the one in the on-shell scheme as Mjj. Therefore, one term in large Higgs boson mass expansion 
is also lost. As a result, if the expansions of |Q are taken, the final result can be given with five coefficients in 
both expansions in the large mass case. The formulae can be found in Appendix The mass difference expansion 
requires an independent calculation of the on-shell propagator diagrams and the result can be found in Appendix |^. 
The numeric results can be found in Fig. ^. It should be stressed, that it was checked that the exact analytic result 
without expansions obtained by the translation procedure described above and by an explicit renormalisation in the 
on-shell scheme are the same. 

It is interesting to consider the transition between the schemes performed purely numerically. In Fig. H, the solid 
curve represents the one-loop correction as well as the sum of the one- and two-loop corrections. The fact that they are 
indiscernible in this scale is due to their relative smallness. The most reliable way of obtaining the correction (apart 
from the exact method) is to take the one-loop result and substitute the MS parameters only in the normalisation 
in Eq. |5^, whereas the masses in (Arj^^j)^^ should be left in the on-shell scheme. This is shown in the curve 2). If 
one, however, simply takes the whole invariant and substitutes all of the MS parameters, then curve 1) is obtained, 
which diverges strongly for Higgs boson masses larger than about 250 GeV. It turns out that the sum of the one- and 
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FIG. 8: Numeric translation of Ar form the MS scheme to the on-shell scheme vs. the exact resuh (solid line). 

two- loop corrections does not reduce substantially the scheme dependence, as shown by curve 3), where the correction 
up to two-loop order in the MS scheme has been given for MS parameters translated from on-shell values using Eqs. |5^ 
to ^. 

VI. COMPUTATIONAL METHODS 

The calculation of the bosonic corrections to the muon life time is a relatively complex task. The number of 
Feynman diagrams to be calculated is around 5000 in gauge. This makes it necessary to use automated software. 



A. Software and checks 



The first step of the calculation is the generation of diagrams. Several systems are presently available. Obviously 
each differs in its easy of use, speed and design concepts. 

The on-shell calculation was based on the C-I--I- library DiaGen |24|. It generates all diagrams together with all 
necessary counterterms. The main advantage of this software is the speed, since all of the diagrams were generated 
in a few seconds, thus making the generation phase a negligible part of the calculation. 

Alternatively, for the calculation with the tadpoles the input generator DIANA |^ has been applied. We note 
that according to the rules, given in Section IV, no counterterm diagrams should be generated. They are all taken 
into account by the multiplicative renormalisation. 

The diagrams to be evaluated can be divided into two broad classes. First are these which can be reduced to 
vacuum bubbles. Here, partial integration identities supplied with analytical formulae can be used. 

The second more complicated problem is the evaluation of the two-loop two-point functions at non-vanishing 
external momentum (at the values = and = My^ in our case). From the several possibilities two different 
algorithms have been used to deal with these diagrams. 

The algorithm described in ||2^ has been chosen because of its simplicity. As an end result of the tensor reduction 
scalar two-loop propagator integrals are obtained. A high precision numerical evaluation of these is currently possible 
with one dimensional integral representations ]2^ . To this end CH — h programs were used based on the library S2LSE 
p8[ . For large scale differences which occur when the Higgs mass is much above the masses of the W and the Z 
boson double precision turns out to be insufficient. An easy way to see it is to remark that the individual terms in 
the result can behave as whereas due to the screening theorem the whole result behaves at most as Mfj. 
For a Higgs boson mass of the order of 1 TeV, this means that cancellations of the order of 10^ will have to occur. 
If we combine this with the fact that in double precision some of the integrals can only be evaluated to 5 digits, the 
numerical instability becomes apparent. A way out of this problem on 32 bit machines is to use software emulated 
quadruple precision. Of course this signifies an important drop in effectiveness. In practice, the software runs about 
20 times slower. Ten times are due to the use of software emulation for arithmetical operations and two to more 
integration points which are needed for higher precision. On present GHz processors, the evaluation of a single point 
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FIG. 9; A choice of formally heavy lines in the large mass expansion of two W boson propagator topologies. 

of the final result requires around 20s and a conservative estimate of the error over the whole range of Higgs boson 
mass from 100 GeV to 1 TeV is four digits. 

Alternatively to the numerical method, we used also the semianalytic method of expansions (see next subsection). 
In this case the huge cancellations mentioned above do not cause any problem. 

The size of the programs written in C++ and in FORM requires stringent checks. A helpful property of the 
bosonic corrections to the propagators is that the value of every single diagram can be obtained rather easily through 
low momentum or large mass expansions. In fact for the Z boson propagators a low momentum expansion up to 
tenth order provided a five digit agreement with the integral representations for each diagram independently and for 
the whole sum. Additionally, we also made an expansion around the point Mh = Mz (see next subsection) and 
got excellent agreement between the numerical and the expanded results. In the case of the W boson propagators 
not all of the diagrams are below threshold. It turns out that 345 contain a photon or a massless ghost line, which 
makes as much as around 160 of them to be either on threshold or infrared divergent. In this case the low momentum 
expansion either fails to converge or converges very slowly. A way out of this is given by large mass expansions. If 
the lines which are to be considered as heavy are chosen in a specific way, then the large mass expansion leads only 
to vacuum bubbles and one-loop propagator diagrams and the convergence is comparable to the case of the Z boson 
propagators. An example choice of the heavy lines for two different topologies is given in Fig. ^. This procedure fails 
only for graphs which represent pure QED corrections to a boson line. In this case however, the result is known 
analytically pi]] . 

Another way of testing the analytical reduction and the diagram generation software is to check the Ward-Takahashi 
identities for the propagators. Here the following relations have been evaluated 

{^zl.L + ^^MzU^zh.) + M'zTI^gIg. - P' {tI^zLY + ^zIl^gIg, = 0, (65) 

p'{<mL-'^Mwn^wG„)+M^n^^l^^.-p^[n^^^^^^^^ = o (66) 

both for on-shell values of the momentum and in an expansion around zero up to third order. Here Gz and Gw stand 
for the neutral and charged would-be Goldstone boson respectively and the subscript "L" denotes longitudinal parts 
of the vector boson self-energies and the scalar vector transitions are given by 

KgM^P^^voAp^), (67) 

where p is the ingoing momentum of the vector boson. 

The combination of the two checks described above, tests the software from the diagram generation to the numerical 
evaluation. An additional test is of course provided by gauge invariance and indeed the calculation was performed 
in the general gauge with three independent gauge parameters. We have observed explicitely the cancellation of 
each of them from the final result and the counterterms. 

Since the bosonic corrections to the propagators in the MS scheme have been evaluated within the large Higgs 
boson mass approach in p3|] a comparison was also possible for the whole result. It turns out that the agreement is 
perfect for Higgs masses running as low as 200 GeV. 

To complete the description of the computational methods, let us note that C++ and FORM were supplied with 
a collection of AWK and Bourne shell scripts managed by several Makefiles. The system prepared in this way runs 
completely automatically from the beginning with diagram generation up to the numerical evaluation with plots. 
Actually, the specificity of the problem allowed to reduce the evaluation time of the whole problem down to only one 
hour and a half, which is rather short for multiloop calculations. 

B. Expansions 

Here we give more details on how the expansions are performed in two different regimes that we considered 
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• in the mass difJerence hz = (M'jj - Mj)/M'^ and 

• in the mass ratio zh — Af| /Mfj . 

The expansion in the mass difference — is especialiy simpie. It is just a Taylor expansion of all Higgs 
propagators and Higgs boson masses in the vertices around Mz- No additional subgraphs are necessary in this case. 
The expansion in the heavy Higgs boson limit is somewhat more involved. It is given by the rules of asymptotic 
expansions Jl3t . 

In addition in the presence of both Mz and Mw we expand in the difference of these masses as well. Indeed 

\, '^ =sin^gv^^ 0.23, (68) 

is a rather small parameter and the convergence of this series is quite fast. This trick was used previously in [ p3[ . 
The advantage of this approach is that in the case of on-shell Green functions all integrals have only one scale. This 
allows one to use the FORM package ONSHELL2 |^ to evaluate these integrals analytically. 

We should also note that to extend the range of the hz expansion we apply the Pade approximation. Throughout this 
paper we use a [3/3] Pade approximant for Ar and [4/4] for the scheme transition formulae. The Pade approximation 
for the Zh series does not work well since this series is nonalternating. 

VII. CONCLUSIONS 

The recent calculation of the two-loop bosonic corrections to Ar performed by two independent groups has been 
described in detail, from the matching onto the Fermi theory to the renormalisation and the explicit results in the 
on-shell and MS schemes. The framework for the evaluation of the Fermi constant Gp based on the low energy 
factorisation theorem has been constructed. It allows one to compute Gf as a Wilson coefficient in a simple manner. 
This approach is general and is also applicable to other low energy quantities. 

A comparison of different expansions and numerical methods has been given. It has been proven that in the wide 
range of Higgs boson masses expansions provide as much precision as needed and cover the whole region of interest. 
The only problematic region, however, is connected to the thresholds for W and Z boson pair production. If the 
Higgs boson was indeed found in this range, then a precise result could also be obtained with expansions but this 
time of the threshold type. The coincidence of the numerical and analytical results serves as a strong check of the 
calculation. 

The accuracy of the numerical transformation between MS and on-schemes has been tested. It is shown that for 
the Higgs boson masses larger than ~ 250 GeV the two loop correction does not reduce the scheme dependence which 
can be explained by huge cancelations of large terms during the transition procedure. 
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APPENDIX A: [/(I) WARD IDENTITY AND THE RENORMALISATION OF CHARGE 

In this appendix we present a derivation of the relation between the charge renormalization constant and different 
wave function renormalization constants valid to all orders of perturbation theory. The derivation is based on the use 
of the U (1) Ward-Takahashi identity for the weak hypercharge gauge group. To begin with, let us take the bare U{1) 
gauge boson field 5° and rewrite it in terms of mass eigenstates 

Bl^c'^^Al + sUl- (Al) 
Here = cos 0^ and = sin 6^, are bare values of cosine and sine of Weinberg angle. 
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In the next step we express our bare gauge boson fields through the renormahzed ones 



(Zf {cwA^ + swZ^} = l^Z^zZ^, + {Z^^f^ A^^^sl. \^{Zzzf'^ Z^ + ^^z^^m} • (A2) 



Now taking the coefficient in front of in the equation above we have 



( 7B\^/'^ _ (7 a/2 . 1 

) Cw — l,^77J + -s^ 



(A3) 



To complete the derivation we need to relate the renormalisation constant to the charge renormalization constant. 
The electric charge is related to the weak hypercharge via the following equation 



e = 9iCw = [Z^ ) gic^Y = Z^e = Z^gicw, 



(A4) 



— 1/2 

where we have made use of U{\) Ward-Takahashi identity Zg-^ = (^^) • Now we can easily deduce, that 



.0 

w 

Cw 



(A5) 



Substituting this relation into Eq. A3 we have 



1 — (-^77) + -^T^Zz^ 



(A6) 



Using this final relation one can considerably simplify the calculation of the on-shell charge renormalization constant 
and avoid dealing with infrared rearrangement while computing the three-point Green function. 

APPENDIX B: LARGE HIGGS BOSON MASS EXPANSION OF Ar^,^' IN THE ON-SHELL SCHEME 

bos 

In this appendix, the on-shcU rcnormalised ArJ^^^ is given in a twofold expansion, in the large Higgs boson mass 
and in the mass difference between the W and the Z boson. The number of terms is consistent with the result 
as explained in Section 0. The leading behaviour both in the Higgs mass and in the sine of the Weinberg angle has 
been factorised out. 
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The occurring transcendental numbers are 



- 7!' 



(B2) 



4 Cl2(7r/3) 
9 V3 



0.2604341376321620989557291432080308... 



while zh — M^/Alfj. Note that the leading term in the Higgs boson mass can be resumed in siii^ 6w to give the 
behaviour 
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The expansion coefficients read (the first four of them expanded to the order 0{z^) were already published in Q) 
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APPENDIX C: MASS DIFFERENCE EXPANSION OF Ar^^j^ IN THE ON-SHELL SCHEME 

(2) 

The correction Ar^^^ in the on-shell scheme for Higgs masses in the vicinity of the Z boson mass is correctly 
described by an expansion in the mass difference between the Higgs boson and the Z boson and in the mass difference 
between the W and Z bosons. The series below contains five terms in both variables 
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The transcendental numbers are the same as in the previous section. The lack of logarithms of mass ratios follows 
from the fact that a Taylor series in the mass difference does not lead to any infrared problems. The variable hz 
denotes {Mfj — M'^)/M^. The first four of the coefhcients were already published in |^ 
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APPENDIX D: LARGE HIGGS BOSON MASS EXPANSION OF Ar^^j^ IN THE MS SCHEME 

In this appendix, Arf^^^ renormalised in the MS scheme is presented as a twofold expansion in the large Higgs boson 
mass and in the mass difference between the W and the Z boson. The expansion is parametrised as follows 



The parameters, i.e. masses and the coupling constant are in the MS scheme. Apart from the numbers Eq. B2, it is 
assumed that ln(m^^) = ln(TO^ /i being the MS renormalisation scale. 
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APPENDIX E: MASS DIFFERENCE EXPANSION OF Ar^^^s IN THE MS SCHEME 

The correction in the MS scheme is given by six coefficients in the double expansion in the mass differences between 
the W and Z bosons and between the Higgs boson and the Z bosons 
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All parameters arc in the MS scheme and hz = {rnjj — m\)/m\. Note also that the logaritms contain the renormal- 
isation scale as \n{m%) = ln(m|/^^). 
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